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Abstract 

We show that a liquid mixture in the thermodynamically stable homogeneous phase can undergo a phase- 
separation transition when rotated at sufficiently high frequency ui. This phase-transition is different from the 
usual case where two liquids are immiscible or where the slow sedimentation process of one component (e.g. a 
polymer) is accelerated due to centrifugation. For a binary mixture, the main coupling is due to a term oc Ap(u;r) 2 , 
where Ap is the difference between the two liquid densities and r the distance from the rotation axis. Below the 
critical temperature there is a critical rotation frequency ui c , below which smooth density gradients occur. When 
uj > uj c , we find a sharp interface between the low density liquid close to the center of the centrifuge and a high 
density liquid far from the center. These findings may be relevant to various separation processes and to the 
control of chemical reactions, in particular their kinetics. 

Separation de phases des liquides miscibles dans une centrifugeuse 
Resume 

Nous prouvons qu'un melange de liquides dans une phase homogene thermodynamiquement stable peut subir 
une transition de phase une fois centrifuge a une frequence uj suffisamment elevee. Cette transition de phase est 
differente du cas habituel ou deux liquides sont non-miscibles ou du cas ou un processus lent de sedimentation d'un 
composant (par exemple un polymere) est accelere par une centrifugation. Pour un melange binaire, le couplagc 
principal est du a un terme oc Ap(u>r) 2 , ou Ap est la difference des densites entre les deux liquides et r est la 
distance de l'axe de rotation. Au-dessous de la temperature critique, il y a une rotation critique de frequence uj c ; 
au-dessous d'uj c , de faibles gradients de densite prennent naissance. Quand uj > uj c , nous trouvons une interface 
mince entre le liquide de faible densite pres du centre de la centrifugeuse et un liquide a haute densite loin du 
centre. Ces resultats peuvent etre appropries a divers processus de separation et au controle de reactions chimiques, 
en particulier de leur cinetique. 
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1. Introduction 



The phase of a liquid mixture typically depends on the ambient temperature, pressure, and composition, but 
may also be influenced by external fields such as magnetic and electric fields [1,2,3,4,5], surface fields [6,7], or 
gravity [8,9,10,11,12,13]. The gravitational field couples to the density of the different components, and thus 
tends to separate them [14,15,16]. Since this coupling is rather weak, gravity is mainly associated with the slow 
sedimentation process of an already immiscible mixture. A common way to accelerate this sedimentation is by the 
use of centrifugation, where the high achievable rotation frequencies permit an effective acceleration highly superior 
to the simple gravitation case [17,18]. Since centrifugation is very effective as an aid to separation of substances 
which are not thermodynamically mixed, it is rarely used on thermodynamically homogeneous mixtures. Indeed, 
as we show below, when a homogeneous mixture is rotating around some symmetry axis, only small composition 
gradients appear in the radial direction with respect to this axis [18]. However, and this is a surprising result of 
our analysis, if the mixture is rotating rapidly enough, above a threshold value, the mixture composition exhibits 
a sharp jump from one composition to another; this sharp change is the signature of a phase-separation. This 
phase separation should be accessible for many mixtures using modern ultracentrifugation techniques. A simple 
theory is developed, analyzed and discussed below. 



2. Theory and Results 

The mixture under consideration here is binary and symmetrical, consisting of two polymers A and B, each 
made up of N identical monomers of volume vq. The case of simple liquids is obtained when N — 1. The local 
volume fraction of the A molecules is <f> (0 < (f> < 1), and the dimensionless Flory parameter \ ~ 1/T characterizes 
the repulsion between the species. The average mixture composition is denoted <f)o and the mass density is po. We 
further assume a linear constitutive relation between composition and density, namely 

p(4>) = po + Ap(0 - 0o ) (1) 

where Ap = pA ~ Pb, and pA and pB are the densities of the A and B components, respectively. Q] The starting 
point of our formulation for a mixture rotating in an angular frequency u) around some axes is the following 
free-energy density 

/ = — / m (0) - ip(^)(cr) 2 (2) 
vo 2 

In the above / m is the Flory-Huggins free energy of mixing, 

/»W = ^01n(0 + ^(1 - 0) m(l - 0) + x4>0- - 0) + ^ 2 (V^) 2 (3) 

where fcs is the Boltzmann constant and T is the temperature. The last term, proportional to (V</!>) 2 , is due 
to an energy penalty for having composition gradients. It gives rise to a characteristic length scale of phase- 
separated domains scaling as ~ (T c — T)~ 1//2 , diverging at the critical point [19]. In the discussion below, we are 
not restricted to mixtures close to the critical point, and therefore this term is subsequently neglected. The phase- 
diagram derived from this mixing energy has a critical point at (<f> c ,Xc) — (1/2, 2/N). The transition (binodal) 
value of x at a given composition <j) is given by 

^ = 2^i ln (r^) w 

We may write to a good approximation T — B/x, where B is a constant with units of K, and in this case the 
critical temperature is T c = BN/2 and the transition temperature is Ti = Bj\t. We are interested in mixtures in 
the homogeneous phase, at temperatures T below T c but above Tt, and off-critical compositions ((/> 7^ 1/2). The 
second term in Eq. (2), proportional to (ur) 2 , is due to the centrifugal forces. Here r is the distance from the 
rotation axes, and is bounded by the smallest and largest radii: Ri < r < R2- 



Our theory does not depend on the linear relation between mass density and composition. Quadratic and higher order dependences 
lead to qualitatively similar conclusions. 



2 




r[m] 



Figure 1. Composition profiles for a binary liquid mixture confined in a tube with inner radius R± — 10 cm and outer radius i?2 — 40 
cm, as obtained from Eqs. (5) and (6) with Xc < X < Xt (T c > X > T±). At low rotation frequencies uj — 10 s -1 , the high-density 
component is pulled toward the larger radius, and the resulting profile exhibits weak gradients. If the rotation frequency is large 
enough, uj — 80 s" 1 > uj c , a sharp interface separating large and small tp- values appears. As the frequency is further increased to 
uj — 300 s — 1 , the difference between coexisting domains becomes larger, and the front moves to smaller values of r. The horizontal 
dashed line is the lowest possible value of the front location -R m i n [Eq. 7)], corresponding to maximal phase-separation. We used 
JV = 200, v = 2 X 10~ 28 m 3 , p = 1000 Kg/m 3 , Ap = 100 Kg/m 3 , cf> = 0.4, T c = 80°C and the temperature is T - T t - 3.5°K. 

In the presence of Earth's gravity alone, such an inertia term is absent and one has to consider buoyancy 
effects only very close to the consolute point. However, as we will show below and in particular for modern 
centrifuges rotating at high speeds, these inertia effects are important and can modify the thermodynamic aspects 
of separation phenomena. 

We minimize the free-energy in Eq. (2) to obtain the following governing equation for the profile <f>(r): 



kgT 
Nv 
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1-0 
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" ~Ap(o>r) 2 



(i = 



(5) 



fi is the Lagrange multiplier (chemical potential) required to conserve the average composition of a system enclosed 
in a volume V: 



(r)d 3 r 



(6) 



Eq. (5) gives an analytical inverted expression for the profile r as a function of (j>. 

In Fig. 1 we display the calculated profile <j>(r) for several angular frequencies uj. Before centrifugation, the 
mixture is initially homogeneous, and thus X < Xt (T > Tt). The centrifugal term in Eq. (5) is spatially-dependent, 
and this means that when the mixture is rotated, (j>(r) has smooth density changes: high-density material is pulled 
toward large radii whereas small density material is pulled toward small r. This behavior is found for \ < Xo, 
irrespective of uj. A different scenario occurs for Xc < X < Xt, which corresponds to a temperature lower than 
the critical temperature but higher than the binodal. At this regime, if lu is small, smooth composition variations 
occur. However, there exist a critical rotation frequency uj c , above which a phase-separation transition takes place. 
If uj > uj c , cf>(r) has a jump from large values at large r to low values at small r. Note that <f>(r) is not constant in 
each of the domains, but instead has a weak gradient, as is seen in Fig. 1. The equilibrium location of the front 
separating the two domains, R, moves to smaller r values as the rotation frequency is further increased above 
uj c . As R decreases, the difference between the compositions in the two domains increases. R has a lower bound 
denoted i? m i n , which corresponds to the maximal possible separation between the two liquids: pure A component 
is found at large r's and pure B is closer to the centrifuge's axes. It follows that 7? m i n obeys 



Rr, 



R\4>a 



(7) 
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Figure 2. Equilibrium front location R separating high- and low-0 domains as a function of angular frequency uj, as obtained from 
Eqs. (5) and (6). R is smaller than R2 and decreases as u increases above u c ~ 12 s _1 ; it tends to R m i n [Eq. (7)] for very large lo. 
The approximate expression for oj c , Eq. (9), gives u> c — 15 s _1 . All parameters and constants are the same as in Fig. 1. 

In order to obtain an analytical expression for the critical frequency ui c , we continue with the following ap- 
proximation. We expand Eq. (5) to linear order in <j> — cj>o, and substitute in Eq. (6). This allows us to find 



— keT f 

A* Jm 
V 



-Ap^iRl+Rf) 



(8) 



where f' m = df m /d<j>. We notice that at the transition frequency and if (j> < 4> c = 1/2, the largest composition, 
found at r = R2, crosses into the unstable part of the phase-diagram. This means that cu c occurs when <j> = <f>t and 
r — R2 are substituted in Eq. (5). The transition composition (j>t at a given Flory parameter \ ( or temperature 
T = B/x) is given by N\ = ln(0t/(l — <f>t)) /(2<j)t — 1), and hence we obtain the following expression for the 
critical demixing frequency 



4fcsT /£( 



v Ap R\ - R\ 



{<t>t - <h) 



(9) 



The notation = d 2 f m /dcf> 2 have been used. A similar expression exists for lo 2 as a function of T — Tt. 

A plot of lo c as a function of the distance from the transition (binodal) temperature T t is shown in Fig. 3. u> c 
increases as the temperature is elevated above Tt. liu> is larger than the plotted uj c , demixing occurs, otherwise the 
mixture displays only weak composition gradients. The numerical value of ui c is quite low, and certainly achievable 
in conventional laboratory equipment - u c = 60 rad/sec is equivalent to 572 rpm. 



3. Conclusion. 



Centrifugation of mixtures, whether of simple liquids or polymers, is shown here to be an effective means to 
induce a phase-separation transition. As is quite intuitive, composition variations occur at any rotation frequency, 
because of the force which appears due to the different densities of the mixture's components [12,13,16,18]. The 
simple mean-field free-energy model we use reveals a surprising phenomenon: for temperature lower than the 
critical temperature T c , there exists a critical rotation frequency ui c - If the rotation is fast enough, lo > uj c , a 
demixing transition occurs where the high density components is pulled toward the outer radius of the centrifuge, 
lower-density components are pulled in the other direction, and a sharp interface appears. At these high rotation 
frequencies, the centripetal force is large enough to cause large composition differences, which destabilize the 
homogeneous mixture in favor of a separated one. A somewhat similar phenomenon was predicted for liquid 
crystals under a gravitational field [8]. 
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Figure 3. Critical rotation frequency uj c for demixing, as obtained from Eq. (9). T is the temperature and Tt is the transition 
(binodal) temperature for demixing without ccntrifugation. We took 0o — 0.2, and all other parameters are as in Fig. 1. 

We note that demixing is also possible with spatially-varying electric fields. The wedge capacitor, made up of 
two flat electrodes tilted with respect to each other, is the analogue of the centrifuge. An electric field couples to 
the different permittivities of the liquid components Ae = ea — £b- The electrostatic energy density is ~ AeE 2 
and in this geometry the electric field varies like r" 1 [4,11]. We thus find that Ap plays the role of the permittivity 
difference Ae, and lot plays the role of an electric field. The critical voltage in the electrostatic case plays the 
same role as the critical frequency ui c in the present case. 

We obtained numerical solution for a closed rotating container with a fixed composition, as well as an approxi- 
mate analytical expression which compares well with the numerical results. The results we find are not specific to 
the mixing free-energy model utilized, Eq. (3); the transition appears in other models which give a qualitatively 
correct phase-diagram in the absence of centrifugation. In addition, the (V<^>) 2 term of Eq. (3) can be accounted 
for. It is expected to decrease the composition gradients and may be important close to the critical point. Numer- 
ical and analytical expressions for the case of an open system, having a connection to a liquid reservoir, and thus 
having an imposed chemical potential, are easier to obtain and are also at hand. However, to experimental setup 
is not straightforward [20], and we chose to leave it out. 

The demixing behavior described above is quite different from the one commonly observed in laboratory ex- 
periments, where the slow sedimentation of a substance (e.g. polymers) in a solvent is simply accelerated by 
centrifugation. Our treaties show that even thermodynamically-stable mixtures of simple and complex fluids can 
be efficiently separated. It seems that this demixing transition can be used in many technological and chemical 
processes. For example, the reaction kinetics of two or more chemical species in a liquid environment can be 
controlled: if the mixture is homogeneous the reaction takes place everywhere, but when it is demixed the reaction 
can only take place at a thin boundary layer, and thus it has slower kinetics. Variations on this theme can also 
lead to accelerated kinetics: if one chemical species is an inhibitor for a reaction, after phase-separation it will 
have less contact with some of the liquid components, and the reaction can take place at a much higher rate. 

Finally, we should stress that we used the simplest mean field model. This approximation can be relevant 
for systems far from critical point and indeed we have shown that even in such situations phase separation can 
be induced for rotation speeds easily attainable by modern centrifuges. The kinetics of phase separation and of 
formation of concentration profile within the sample seems to be an interesting open question that we leave until 
some experiments are performed. 
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